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A Sunyaev-Zel'dovich take on cluster radio haloes - 1. Global 
scaling and bi-modality using Planck data 
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ABSTRACT 

Giant radio haloes in galaxy clusters are the primary evidence for the existence of relativis- 
tic particles (cosmic rays) and magnetic fields over Mpc scales. Observational tests for the 
different theoretical models explaining their powering mechanism have so far been obtained 
through X-ray selection of clusters, e.g. by comparing cluster X-ray luminosities with ra- 
dio halo power Here we present the first global scaling relations between radio halo power 
and integrated Sunyaev-Zel'dovich (SZ) effect measurements, using the Planck all-sky cluster 
catalog and published radio data. The coiTelation agrees well with previous scaling measure- 
ments based on X-ray data, and offers a more direct probe into the mass dependence inside 
radio haloes. However, we find no strong indication for a bi-modal cluster population split be- 
tween radio halo and radio quiet objects. We discuss the possible causes for this apparent lack 
of bi-modality, and compare the observed slope of the radio-SZ correlation with competing 
theoretical models of radio halo origin. 
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1 INTRODUCTION 

The bulk of the baryonic mass in galaxy clusters exists in the form 
of a low density, diffused ionized gas filling up the space between 
cluster galaxies. This hot intra-cluster medium (ICM, T ~ 2 - 10 
keV) emits in the X-rays through thermal Bremsstrahlung emis- 
sion, and is also observable in the millimeter/sub-millimeter wave- 
lengths through the Sunyaev-Zel'dovich (SZ) effect, which is a 
distortion in the intensity of the Cosmic Microwave Background 
(CMB) radiation caused by the same thermal component (Sunyaev 
& Zel'dovich 1980). Together, these two observables are central to 
the use of galaxy clusters as cosmological probes. 

The ICM is also host to a large population of ultra-relativistic 
particles (cosmic rays) and magnetic fields, seen primarily through 
radio observations. The most spectacular evidence for this non- 
thermal population comes from observations of giant radio haloes, 
which are diffuse sources of radio synchrotron emission extending 
over ~ 1 Mpc scales. The haloes are not associated to any par- 
ticular cluster galaxy, and are morphologically distinct from radio 
mini-haloes (residing in cluster cool cores), radio relics (formed at 
the edge of a merger shock) and radio lobes associated with active 
galactic nuclei. The similarity of their morphology with the ICM 
suggests a correspondence between their powering mechanism and 
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the total cluster mass (e.g. Liang et al. 2000). They are relatively 
rare and are found mostly in clusters showing evidence of ongoing 
mergers. As such, they can prove to be essential in understanding 
cluster merger dynamics and associated heating processes in the 
ICM (see e.g. review by Ferrari et al. 2008). 

Despite their importance, the powering mechanism of these 
giant radio haloes remains uncertain. There are two models for par- 
ticle acceleration in a radio halo volume: the hadronic model which 
uses collisions between cosmic-ray protons and thermal protons for 
generating relativistic electrons (Dennison 1980), and the turbu- 
lence models where the electrons are re-accelerated through MHD 
turbulence in the ICM caused by cluster mergers (Brunetti et al. 
2001, Petrosian 2001). The distinction between these two models 
is partly based upon the observed scaling between radio and X-ray 
power (the latter indicating the total cluster mass), and the fact that 
X-ray selection seems to indicate two distinct populations of clus- 
ters: the radio halo and "radio quiet" ones (e.g Brunetti et al. 2007). 
However, recent discoveries of radio haloes in clusters with very 
low X-ray luminosities (Giovannini et al. 201 1), and the lack of ra- 
dio haloes in some mergers (e.g. Russell et al. 201 1) show that the 
X-ray selection may not be as clean as expected. These new obser- 
vations and the underlying large scatter in the Lx - /"radio correlation 
suggest that a new observational window on the selection and mass 
estimation of clusters harboring radio haloes can bring some much 
needed clarity. 
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One further reason for expecting a robust correlation between 
radio power and SZ is the timescale argument: the boost in the 
X-ray luminosity during mergers happens in a relatively short 
timescale, compared to the gas thermalization in a modified poten- 
tial well producing a more gradual and moderate increase in the SZ 
signal (Poole et al. 2007, Wik et al. 2008). This should correspond 
better with the radio halo time scale (~ 1 Gyr), derived from the 
spatial extent of the haloes. The integrated SZ signal is also a more 
robust indicator of cluster mass than the X-ray luminosity, irrespec- 
tive of cluster dynamical state (e.g. Motl et al. 2005, Nagai 2006). 
Thus SZ-selection might be able to find radio haloes in late mergers 
and other massive systems which are left out in X-ray selection. 

This letter presents the first radio-SZ correlation for clusters 
with radio haloes. The radio data is a compilation of published 
results, and the SZ measurements are taken from the Planck all- 
sky cluster catalog (Planck collaboration 201 1). All results are de- 
rived using the ACDM concordance cosmology with Q.^ = 0.26, 
Sl\ = 0.74 and Ha = 1\ km s"' Mpc"'. The quantity Fsz is used 
throughout to denote the intrinsic Compton F-parameter for a clus- 
ter, Yd\, where dA is its angular diameter distance. 

2 RADIO & SZ DATA SETS 

We do not attempt to define a new comprehensive sample for this 
work, rather use a set of available cluster catalogs with radio halo 
detections and non-detections to probe the robustness of the radio- 
SZ scaling. Published radio error estimates often ignore system- 
atic effects like flux loss in interferometric imaging and contribu- 
tion from unresolved point sources, which in turn create an over- 
estimation of the intrinsic scatter. Since the present work is mainly 
concerned with the mean slope of the radio-SZ scaling and not its 
dispersion, error underestimation in the literature will not affect the 
results as long as the measurements are unbiased. 

The radio catalogs can be divided into two groups: those with 
and without a listing of non-detections. A comprehensive sample in 
the former category is by Giovannini et al. (2009, hereafter G09), 
presenting results at 1.4 GHz for z < 0.4 clusters. Potentially prob- 
lematic are its mixing of radio halos and mini halos, and not sep- 
arating contributions from radio relics. More critically, the sizes 
of the radio haloes are approximated by the observed largest lin- 
ear sizes (LLS), which is not a good approximation for radio halo 
diameter. To address the latter issue we use a smaller subsample 
by Cassano et al. (2007, hereafter C07), which provides a better 
measurement of radio halo sizes by averaging their minimum and 
maximum extensions. 

The most systematic study of radio halo non-detections in an 
X-ray selected sample is by Venturi et al. (2008), using GMRT ob- 
servation at 610 GHz. However, this sample contains too few clus- 
ters which have Planck SZ measurements to obtain any robust cor- 
relation. We therefore use the compilation by Brunetti et al. (2009, 
hereafter B09), which lists GMRT results scaled to 1.4 GHz with 
other unambiguous halo detections. The non-detection upper lim- 
its were obtained by simulating fake radio halos in the GMRT data 
and scaling to 1.4 GHz by using a = 1.3, a typical spectral index 
for radio haloes. Our final sample is from Rudnick & Lemmerman 
(2009, hereafter R09) who re-analyzed WENSS survey data at 327 
MHz for an X-ray selected sample. The shallowness of WENSS 
data makes R09 ineffective in testing bi-modality, as the 3(T upper 
limits ai'e not sufficiently below the detection level. Its use is mainly 



Table 1. Regression coefficients for the scaling relation log(Py) = A + 
B log(ysz)- The term global implies correlation with the total SZ signal, 
as opposed to that scaled inside the halo radius. 



Sample 


sub-categoiy 


B (slope) 


A (norm.) 


G09 


global 


1.84 ±0.38 


31.3 ± 1.4 




inside LLS 


0.95 ±0.14 


28.8 ± 0.5 


C07 


global 


1.88 ±0.24 


31.4 ±0.8 




inside Rh 


1.17±0.18 


29.7 ± 0.8 


B09 


global, haloes only 


2.03 ± 0.28 


32.1 ± 1.0 




+ non-detections 


2.41 ± 0.44 


33.4 ± 1.6 


R09 


global, haloes only 


0.81 ±0.36 


28.1 ± 1.4 




+ non-detections 


1.38 ±0.43 


29.8 ± 1.8 



Samples: G09=Giovannini et al. 2009; C07=Cassano et al. 2007; 
B09=Bmnetti et al. 2009; R09=Rudnick & Lemmerman 2009. 



limited to testing possible changes in the scaling law at lower fre- 
quencies. 

The Sunyaev-Zel'dovich effect measurements are taken from 
the Planck ESZ catalog (Planck collaboration 2011). This all-sky 
cluster catalog provides a list of 189 objects in the highest mass 
range out to a redshift z ~ 0.6, selected al S /N > 6 from the first 
year survey data. Out of these, 22 are either new cluster candidates 
or have no X-ray data. The remaining 167 clusters, spanning a red- 
shift range 0.01 < z < 0.55, are cross-correlated against the radio 
catalogs. All radio halo clusters therefore have an R^m estimate in 
the Planck catalog obtained from the Lx - M500 relation. This is 
used to model the pressure profile for each cluster when scaling 
their integrated SZ signal, Ksz. between different radii. 

Proper regression analysis between the radio and SZ observ- 
ables is fundamental to this work. We must allow for measurement 
errors and intrinsic scatter in both observables, which makes the re- 
gression analysis non-trivial (see Kelly 2007). We use the publicly 
available idl code by Kelly to perform the regression analysis us- 
ing a Bayesian approach. An important advantages of this method 
is the provision for including non-detections. 



3 RESULTS 

3.1 Radio-SZ scaling 

The radio-SZ scaling relation is obtained by performing linear re- 
gression in log-space: log(Py) = A + 61og(ysz). The normalization 
A and slope B are obtained from the Markov Chains, as well as the 
intrinsic scatter crioi,/>|ioi,>-. A summary of the results are given in 
Table 1 . The first correlation example is from the G09 sample, com- 
paring the radio power with the global SZ signal (Pig.\l\left). Out 
of 32 objects in this sample 24 have Planck counterparts. The mean 
slope is 1 .84 ± 0.38. There is a lot of scatter in this correlation, with 
mean scatter 0.45 dex, i.e. roughly a factor ~ 2.8. Much of this scat- 
ter is driven by the low-redshift objects which are under-luminous 
(e.g. A401 and A754). This potentially indicates a systematic bias 
in their total flux and size measurements with interferometers. Only 
A1351 stands out as overtly radio luminous for its mass, although 
later revisions of its radio power (Giacintucci et al. 2009) moves it 
closer to the mean value. 

The Planck catalog provides the integrated Y parameter within 
radius SRsoo, obtained from a matched filtering algorithm assuming 
a universal gas pressure profile (see Planck collaboration 201 1). At 
this radius, Y^^' x Y^!!' . This is nearly 3 times the cluster virial 
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Figure 1. Radio-SZ correlation for Planck detected clusters with radio haloes. Left — Result from the G09 sample, correlating radio halo power against the 
total Fsz signal inside 5i?5oo. Filled symbols correspond to clusters at z > 0.2, and open symbols are for lower redshifts. The shaded area marks the 2a" or 
95% confidence region (some name labels are omitted for claiity). Right - The same G09 clusters after scahng the total SZ signal to inside the haloes' largest 
linear dimensions, resulting in a much flatter slope with reduced scatter. (Inset) Result from the C07 sample, with scaled SZ signal inside their quoted radio 
halo radius, Rfj. Mean slope is 1.17 ± 0.18, which is used for comparison with theoretical models. 



radius, and much larger than the extent of the radio emitting re- 
gions. Therefore, a tighter correlation can be expected if this global 
SZ signal is scaled down to that inside the radio halo volume. We 
do this conversion by assuming the universal pressure profile of Ar- 
naud et al. (2010), as also used by the Planck team. In particular, 
the best fit profile for mergers/disturbed clusters from the appendix 
of Arnaud et al. (2010) is used, but the difference is negligible if 
the mean profile is used instead. This scaling of the SZ signal in- 
side LLS changes the results significantly. Correlation between ra- 
dio and SZ powers inside the halo volume becomes consistent with 
a linear relation, with mean slope 0.95 ± 0.14 (Fig\l\right) and a 
reduced mean intrinsic scatter in radio (0.35 dex). 

The largest linear sizes are in general not a good approxima- 
tion for radio halo diameters, so the above analysis is repeated with 
the CO? sample using their revised measurement of halo radius, Rh- 
The slope for the global correlation with this sample is 1.88 ± 0.24, 
whereas after scaling the SZ signal inside Rh it becomes 1.17±0.18, 
with mean intrinsic scatter 0.28 dex (Fig[T] inset). Although this 
is statistically fully consistent with the scaled result of the G09 
sample inside the LLS, we use this slightly super-linear correla- 
tion when making comparisons with theoretical models, due to the 
better definition of halo radius. We emphasize that from the current 
analysis using radio data from the literature, a linear correspon- 
dence between radio and SZ power is a fully valid result. 

The R09 sample at 327 MHz indicates a flattening of the cor- 
relation slope at lower frequencies: the best fit value is 0.81 ± 0.36 
when considering the halo sample, with a scatter of only 0.21 dex. 
The large flux uncertainties (and correspondingly low intrinsic scat- 
ter) reflect on the shallowness of the WENSS data, which is more 
than an order of magnitude less sensitive compared to typical VLA 
measurements scaled to its frequency. However, the method used 
by R09 to detect haloes (and place upper limits) based on simulat- 
ing sources in control regions safeguards against flux underestima- 
tion bias, which will otherwise occur in a visual inspection. There 
can be a residual bias from fluxes associated with small scale struc- 
tures that are not recovered. If non-detection upper limits in R09 
are included in the correlation, then the slope becomes 1.38 ± 0.43, 
which is consistent at Icr with the scaling result at 1.4 GHz. 



3.2 Lack of strong bi-modality 

To test whether there are two distinct populations of clusters: those 
hosting powerful radio haloes and those without, we use the B09 
sample. Actual flux measurements at 1.4 GHz are used for the 
clusters A209 and RXCJ1314 (Giovannini et al. 2009), instead of 
the extrapolated values from 610 MHz as given by B09. All non- 
detection upper limits are extrapolation of simulation results at 610 
MHz using halo spectral index a = 1.3. Regression analysis for 
the two cases (with and without halo non-detections) yields slopes 
which are statistically consistent with each other, even though a bi- 
modal division appears to be emerging (Pig.^left). Significantly, 
we do not find high-Ksz objects with radio non-detections, as in 
the case of highly X-ray luminous "radio quiet" cool core clusters. 
But the small number of non-detections makes it difficult to con- 
clude whether the bi-modality is weaker or non-existent. All non- 
detection upper limits lie below the 95% confidence interval of the 
halo-only correlation, suggesting clusters with upper limits are gen- 
erally radio under-luminous. This can be partly redshift-driven, as 
the samples are not SZ-complete. 

An alternative, although not independent, way to visualize this 
result is to correlate the radio power directly against cluster gas 
mass. Rough estimates of Mg^s inside R^qq were obtained by divid- 
ing the Planck Ksz values by the mean X-ray temperatures taken 
from the literature. The global Ksz measurements of Planck are 
scaled to that inside Rjog using the universal pressure profile, as 
this is the radius within which X-ray temperatures are typically ob- 
tained. The result is similar to the Pi 4 - Ygz correlation, lacking 
a strong bi-modal division (Fig[2l right panel). The scatter is in- 
creased by roughly 30% (from 0.5 dex to 0.6 dex), in line with the 
expectation that Fsz is a lower scatter mass proxy. The four non- 
detection clusters have generally lower Tx values as compared to 
the halo detections (median Tx is 7.3 keV compared to 8.7 keV). It 
should be made clear that by taking Tx estimates from literature we 
ignore potential errors due to non-uniform radius for extracting Tx, 
systematic differences between XMM-Newton and Chandra mea- 
surements, etc. However, the mean slope for the P1.4 - Mga,, scaling 
relation, 3.2 ± 0.7 with the full sample, is consistent with the global 
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Figure 2. Test of bi-modality with non-detections from tlie B09 sample. Left - Correlation of radio power against the total SZ signal: filled symbols are 
for z > 0.2 cluster haloes, and open symbols are for lower redshifts. All non-detections are at ; > 0.2 and are extrapolated from 610 MHz data. The only 
mini-halo (A2390) in the sample is marked by the orange square. The short-dashed line con'esponds to the fit for haloes only, and the long-dashed line when 
non-detections are included. Filled regions mark the 95% confidence intervals. Right — Congelation of radio halo power with the gas mass. Symbols and lines 
have the same meaning as in the left panel. 



mass scaling derived from the Y - M relation (see 34. It . indicating 
that no additional biases are incurred while using this non-uniform 
selection of X-ray temperatures. 

A tentative argument for a selection bias in X-ray complete 
samples and the ensuing bi-modality can be given by comparing 
the relative frequency with which radio haloes and non-detections 
occur in the Planck catalog. In Venturi et al. (2008), GMRT data 
were obtained for a complete X-ray selected sample, with 6 de- 
tection of radio haloes plus 20 non-detections. The Planck catalog 
contains 5 out of these 6 halo clusters, but only 4 out of 20 non- 
detection clusters. For the B09 sample this ratio is 16 out of 21 
radio halo clusters and 4 out of 20 non-detections (the same non- 
detections as in the Venturi et al. sample). Since the Planck catalog 
should not have a significant bias towards mergers, this provides 
an indirect evidence for our hypothesis that being SZ-bright (hence 
massive) is a better indicator for clusters hosting radio haloes, as 
opposed to being X-ray luminous. Even though the R09 sample is 
too shallow to directly test bi-modality, it interestingly follows this 
same trend: Planck reports measurement of 12 out of 14 counter- 
parts for haloes and other diffuse emissions, as opposed to only 15 
out of 58 counterparts for non-detections. 



4 THEORETICAL CONSIDERATIONS 

4.1 Mass scaling of the radio halo power 

The independent variable, isz, is defined as the integral of the total 
pressure in a spherical volume, and hence is proportional to the 
total gas mass: 

Ksz = Ydl^ J n,T, dV oc M,,,T, = f,,Mo,T,. (1) 

Here Tg is the mean gas temperature within the integration radius, 
and /gas is the gas-to-mass ratio. Assuming hydrostatic equilib- 
rium and isothermality, the temperature scales to the total mass 
as T, cc Mll^Eiz)-^^ (e.g. Bryan & Norman 1998), where E{z) is 
the ratio of the Hubble parameter at redshift z to its present value. 
Therefore, the scaling between the SZ observable and total mass is 



i'sz£^(z) oc f^^^M^lf . Numerical simulations, analytical models 
and SZ observations indicate that this mass scaling is extremely ro- 
bust, with little scatter over a large range of cluster mass, dynamical 
state or other details of cluster physics (e.g. Motl et al. 2005, Reid 
& Spergel 2006, Andersson et al. 201 1). We thus assume this scal- 
ing to be valid also inside a cluster at different radii, provided that 
the radius is sufficiently large to exclude complex physics at clus- 
ters cores. The large halo sizes measured by C07 {Rn ~ 600 kpc, 
typically of the same order as Risoa), ensures that they encompass 
a representative cluster volume. The E{zT'^^ factor for self-similar 
evolution changes the scaling results only marginally, well within 
the statistical errors. The gas mass fraction, /„as, has a weak depen- 
dence on cluster mass: /„as cx M^^j^* (Bonamente et al. 2008, Sun et 
al. 2009). Assuming the same mass dependence of for all radii, 
we therefore obtain 

P,, cc M^^^^\ (2) 

In the above, Mh is the total mass inside radio haloes, and Myir 
is the cluster virial mass which scales linearly with M,o,(< SRsoo)- 
The scaling index inside haloes is in good agreement with previous 
X-ray hydrostatic mass estimates (e.g. Cassano et al. 2007). The 
global scaling with total cluster mass can be a useful parameter for 
estimating radio halo statistics, particularly in simulations. 

The radio halo sizes are known to scale non-linearly with 
cluster radius, in a break from self-similarity (Kempner & Sarazin 
2001, Cassano et al. 2007). Indeed, using the X-ray derived R<,qo 
measurements from the Planck catalog, we obtain the empirical re- 
lation Rh <x ^500*""^ with the C07 sample, consistent with the esti- 
mate by C07 using Lx-My,^ scaling relation (R/j oc R^-^-''-^). A con- 
sequence of this rapid increase in radius is a drop of the mean gas 
density inside haloes with increasing halo mass. Our observed scal- 
ing between the halo radius and scaled SZ signal, R^ oc Y^-^'--^'"^ , 
implies that the mean gas density («;/) scales down as roughly 
riH 0^ r^"', or assuming thermal equilibrium inside haloes, as 
hu oc M^"-^. This brings the observed non self-similar scaling be- 
tween Rh and /fvir in conformity with the mass scaling in Eq.l|2j. It 
is worth mentioning at this point that radio halo size measurements 
with insufficient S/N will tend to show a steeper dependence on lu- 
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minosity than the true scahng, since only the bright central regions 
will be picked up. 

4.2 Comparison with radio/X-ray scaling 

There is some confusion in the literature about the exact power 
in the X-ray/radio scaling: reported values using the luminosity in 
the soft X-ray band range between Py oc -L^[Q"f_2 4j (Brunetti et al. 
2007, Kushnir et al. 2009). Using the regression method used in 
this work we find the scaling index in the middle of this range, 
e.g. from the B09 sample the mean slope for log Py — log Lx[aA-2A] 
correlation is 1.80 ± 0.21, with mean intrinsic scatter 0.3 dex. 
The mass-luminosity relation for the X-ray soft band is well- 
established observationally. We use the result given by Zhang et al. 
(201 1) for disturbed clusters in the HIFLUGCS sample: La:[o.5-2] oc 
[^Bas./i5oo^(')] ' where the luminosities are core corrected. 

This combined with the weak mass dependence of /gas produces 
a mass scaling of radio power as /'1.4 oc M^^, which is much shal- 
lower than the virial mass scaling obtained in Eq.Q but roughly 
consistent with the halo mass power law. This indicates that the 
global X-ray emission acts as a relatively good proxy for radio halo 
masses due to its peaked profile, as most of the X-ray flux comes 
from within a radius that is < R^. 

4.3 Expectations from theoretical models 

The hadronic model for radio synchrotron emission postulates that 
electrons at ultra-relativistic energies are produced in the ICM by 
p-p collisions between cosmic ray protons and thermal protons (see 
review by Ensslin et al. 201 1 and references therein). For estimating 
the scaling relation between radio halo power and cluster mass, we 
follow the formulation by Kushnir et al. (2009). In this model, the 
total radio power is the volume integral of the cosmic ray energy 
density (ecR = ^cr n kTg) and the hadronic interaction rate (Tpp ~ 
n o-pp): 

= J Tpp' ecR dV ~ XcR n- kT, o-pp fs R^. (3) 

Here Xcr is the ratio between cosmic ray pressure and thermal pres- 
sure, n is the gas density, (Xpp is the p-p collision cross-section, 
is the volume filling factor for magnetic fields, and R^ is the halo 
radius. In the second step we have assumed the magnetic field en- 
ergy density to be much larger than the CMB energy density, B » 
BcMB ~ 3.2(1 + zfpG. Considering that Ysz{< Rh) ~ n kT, Rl, we 
thus obtain: Py/Ysz ^cr fs n CTpp. Therefore, if the cosmic ray 
fraction and mean density do not depend on the halo mass, we re- 
cover the observed linear dependence between radio power and the 
SZ signal inside haloes. The latter assumption, however, is in con- 
flict with our observed scaling of mean gas density 04.lt . which 
actually drops with increasing halo mass. Another potential prob- 
lem is the assumption of strong magnetic fields, Bff » Scmb> over 
the entire halo volume. A likely scenario with the hadronic model 
would therefore be to assume a more clumpy radio emission, where 
the regions contributing most of the radio power have constant den- 
sities and strong magnetic fields. 

In the turbulent re-acceleration model a pre-existing popula- 
tion of electrons at lower energies gets re-accelerated by merger in- 
duced turbulence (see review by Ferrari et al. 2008 and references 
therein). The powering mechanism of radio haloes is complex, but 
for a simple estimate we can follow the formulation by Cassano 



& Brunetti (2005) and Cassano et al. (2007). In their model, the 
energy injection rate from turbulence depends on the mean density 
and velocity dispersion inside the radio haloes; e, oc n cr^, where 
cr^ = GMh /Rh- The total power of a radio halo is then 

r Mh (rl 

Py~ \ s, (r,„/r,h) dv oc " (4) 

where F is the turbulence damping rate transferring energy to the 
particles, and the function Tiz, M^, B^) is constant in the asymp- 
totic limit of strong magnetic fields. Thus to the first approxima- 
tion, Eq.l[4ll implies P,. oc YhTI^~, in good agreement with the 
slightly super-linear slope inside haloes seen from the C07 sam- 
ple (/'1.4 oc y^' "±'' '**). However, using the definition of cth and the 
observed scaling between halo mass and radius, we find a mass de- 
pendence slightly shallower than observed: Py oc Mjf^'--^, which is 
still consistent with Eq.l|2j. This may indicate a preference for more 
realistic field strengths, e.g. figure 2 in C07 suggests approximately 
!F oc M^"^ if the mean field strength is of the order 5-6 /jG in- 
side a radio halo of mass M/j ~ 10'*^ Mq, assuming B/j oc M^^. A 
shallower Bh - Mh relation will correspondingly imply an weaker 
field to explain the observed P]^- scaling. 

5 CONCLUSIONS 

In this letter we presented the first radio-SZ correlation results 
for clusters hosting radio haloes, using published radio data and 
the Planck SZ catalog. There is a clear correspondence between 
these two thermal and non-thermal components as expected from 
the well-established radio/X-ray correlation. On the other hand, we 
found no strong bi-modal division in the cluster population split be- 
tween radio halo and "radio quiet" objects. The halo non-detection 
clusters are generally radio under-luminous, but their occurrence in 
the Planck catalog is much less frequent as compared to all X-ray 
complete samples, and as such we can not conclude whether the 
bi-modality is weaker or non-existent when measured against SZ. 
A likely explanation for this difference can be that the bi-modality 
seen in the Lx selection comes from a bias towards lower mass cool 
core systems (which are radio quiet), whereas SZ selection picks up 
the most massive systems irrespective of their dynamical states. A 
forthcoming work will purport to test this hypothesis using a com- 
plete SZ-selected sample. 

The radio-SZ correlation results were compared with the 
simplified theoretical predictions from hadronic and turbulent re- 
acceleration models. Even though the observed correlation power 
can be explained from both these models under certain assump- 
tions, the turbulent re-acceleration model can be considered abetter 
fit to the data given the simple formulations used in this letter. The 
difference between the global radio-SZ scaling and the one within 
the halo volume is explainable from the non-linear scaling between 
radio halo mass and total cluster mass. An indicative flattening of 
the correlation slope was observed when considering a cluster sam- 
ple at 327 MHz, but the result became consistent with 1.4 GHz 
observations when the haloes and non-detections were considered 
together as one population. 
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